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Abstract 

Observations suggest, that there may be periods in the history of the universe, including the present 
one, in which its evolution is driven by scalar fields. This paper is concerned with the solution of the 
c;vohition equations for a spatially flat universe driven by a single scalar field. Some general theorems 
relevant to the cosmology of these models are presented, and several approaches to solve the equations 
are discussed. For some potentials special exact solutions can be found, and for the case of exponential 
potentials the complete solution is rederived in a new parametrization. For the general case solutions 
arc constructed in terms of a power series expansion in the field. The issue of double-valuedness of 
such a series expansion in case of oscillating fields with turning points is addressed and resolved. 



1 Cosmic scalar fields 



It is widely accepted that scalar fields can drive the cosmic expansion. More in particular, 
scalar fields can possibly account both for an early period of inflation to explain the large-scale 
homogeneity and isotropy of the universe [H |2] , and for the observed accelerated expansion of 
the universe in more recent times [Sj S] . One of the most widely studied scenarios for inflation 
is the minimal chaotic inflation model [5] , in which a scalar field moving in a potential creates 
a dynamical form of dark energy that makes the universe expand. In a similar way cosmic 
scalar fields can be used to model dark energy [6l [7] which drives the accelerated expansion 
of the universe deduced from supernova observations. In view of such potential applications 
the dynamics of scalar fields in a cosmological context is a physically relevant subject of 
investigation O |9] . 

In this paper we investigate simple cosmological models in which a single dynamical scalar 
is minimally coupled to gravity, and the gravitational field is taken to be of the Friedmann- 
Lemaitre-Robertson- Walker type. More in particular, motivated by the observations of the 
cosmic microwave background we take the spatial part of the metric to be flat. Thus the line 
element describing the proper-time for a comoving particle in the universe i^ 

- dr^ = g^^ dx^'dx" = -N^{t) dt^ + a^{t) dx^. (1) 

Here a{t) is the scale factor, whilst N{t) is the lapse function in the ADM formulation of 
General Relativity |10j, allowing one to keep local time reparametrizations as an invariance 
in the description of space-time geometry. Usually one chooses t to be cosmic time such that 
N{t) = 1. This will also be our prefered choice. However, we find it useful to keep N{t) free 
in the derivation of the relevant field equations for reasons to become clear soon. 

The action for a real scalar field if minimally coupled to gravity in General Relativity is 

S = ld^xV^(^-^R-^g>^''d,^d,^-V[^]y (2) 

where R is the Riemann scalar and V[ip] is the scalar potential. Taking the metric of the 
form ([T]) and the scalar field to be spatially constant in this frame: V(/3 = 0, ip{x'^) = '^{t), 
the effective action per unit of co-ordinate volume takes the form 

dt{-^ac? + ^^^-a^NV[^]^. (3) 

Variation of this action provides the relevant equations of motion for the evolution of this 
isotropic and homogeneous model universe. First, the dynamical equation for the scale factor 
is 

I 511 _ 2 d f I da\ f I daV 1/1 dipV v\ ] - 

Next, the dynamical equation for the scalar field is 

1 5S 1 d fld^\^/3 da\ld^^^,^^^^^^ 



a^N 5ip N dt\N dt J \aN dt J N dt 



^In most of this paper we employ Planck units such that c — h — SttG — 1. 
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where the prime denotes a derivative w.r.t. the field 99. Finally, there is a constraint imposed 
by variation of the arbitrary lapse function A^: 

1 5s /I day 1 fid^y 

In this form, both the effective action and the equations of motion are manifestly invariant 
under time reparametrizations t — )■ t' , with N{t) transforming as 

N'{t')dt' = N{t)dt, (7) 

whilst a{t) and ip{t) behave as scalars: 

a'it') = a{t), ^'{t') = ^{t). (8) 

This reparametrization invariance allows us to choose a gauge = 1, such that the equations 
of motion become 

27? + 3i/2 + ^ 02 - F = 0, 

(p + 3Hip + y = 0, (9) 
-3^2 + 1 ^ y = 0, 

where H = d/a is the usual Hubble parameter. It is well-known that these equations are 
redundant to the extent that for i? 7^ the last two equations imply the first one; indeed, for 
H ^ differentiation of the last equation leads -upon use of the middle one- to 

6HH + SHip"^ = ^ 2H = -ip^. (10) 

Adding this to the last equation ([9]) gives back the first one. However, when H = the 
proof fails, and the first equation has to be considered as a separate condition. Actually, the 
last equation ^ obtained from the variation of N is just the reduced hamiltonian constraint 
of general relativity in FLRW space-times (also known as the Wheeler-DeWitt equation) 
which restricts the set of allowed solutions of the equations of motion to those for which the 
hamiltonian vanishes. 

In this context we observe, that the effective action ([3| can be written in a more familiar 
form by defining new dynamical variables {X^,X^) and a time parameter r(i) by 

X^ = V6lna, X' = ip, ^ = \ = e~Vy~^^' . (11) 

dt a-^ 



(12) 



In terms of these variables the effective the effective action becomes [12 

1 /dX^V 1 /dX^^^ 



where 

U[X°, X^] = e^^V[Xi] = a^V[^]. (13) 

Thus the cosmological model ([s]) is mathematically equivalent to that of a relativistic particle 
in Minkowski space moving in a time-dependent scalar potential U. The action (12) is a 
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convenient starting point for a canonical (and quantum) treatment of mini-superspace cos- 
mofogy. In this formulation the hamiltonian constraint in the gauge = 1 takes the simple 
form 

'''"'°^\U'^Umx',x']-». (14) 



2 \ cIt J 2\ dr ^ 

2 Dynamics: general considerations 

Returning to the classical cosmology model described by eqs. ([o]), we observe that the hamil- 
tonian constraint represents a first integral of motion for the system (a, (^), but one which 
can not take arbitrary values: being a first-class constraint the right-hand side of the third 
equation ^ necessarily vanishes, even though the other two equations would be consistent 
with any constant value E such that 

(^-3H^ + ^ip'^ + = E. (15) 

The constraint -E = is the result of local time-reparametrization invariance imposed by 
the gauge variable A^. It follows that one can not impose arbitrary initial conditions for the 
variables (a(t),(^(t)) and their velocities {a(t),ip{t)): any set of initial values is constrained 
hy E = 0. Keeping this in mind, the complete classical dynamics can be derived from the 
hamiltonian constraint and the Klein-Gordon equation for ip, i.e. the second eq. ([9]). 

It is clear that the middle term of the KG equation linear in the velocity ip represents the 
transfer of energy from the scalar field to the scale factor, or vice versa. Note however, that 
it does not imply the breaking of time reversal invariance: under time reversal t — )• — t both 
if and H change sign, with the effect that the equation itself is time-reversal invariant. As a 
result the energy density E/a^ remains constant (and vanishes), and there is no dissipation 
of energy for the combined system of scalar and gravitational degrees of freedom as such. To 
include dissipation by e.g. particle creation, the equations (|9| would have to be modified [H]. 

A related conclusion is, that consistent non-degenerate evolution of the system does not 
necessarily require the scalar potential V to be bounded below: there is already a negative 
contribution to the hamiltonian from gravity, witness the term —3H^ in eq. (15), but the 
constraint E = serves to stabilize the system. 

As observed earlier the dynamics of gravity, as described by the first equation ([o]), only 
follows from the other two equations if H ^ 0; therefore in solving for a and ip from the KG 
equation and the constraint we always have to consider the case H = corresponding to fiat 
space-time separately. The case of Minkowski space-time H = H = is quite straightforward: 
the equations ^ reduce to 

^p^-V=^p^ + V = 0, p + V' = 0. (16) 
It follows that the kinetic and potential energy have to vanish separately: 

p = V = 0, (17) 

and as a result also 

ip = V' = 0. (18) 
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This is possible only if the potential has a stationary point which is also a zero: V = V' = 0. 
In most potentials this will not apply. Note, that of course the evolution of the universe can 
pass through a flat point, where H = but H ^ 0. At such a point 

H = V = -^^' <0, (19) 

i.e. the universal expansion goes through a maximum there. 

In all other cases {H 7^ 0), it suffices to consider the two equations 

ip + 3H^ + V' = 0, ^vj^ + y = 3i/2. 

Now assuming a solution to exist, this solution must at least for finite stretches of time define 
a one-to-one map from t to During such a period one can write |12l 



H{t) = H[ip{t)], H = H'^, (20) 

where the overdot denotes a derivative w.r.t. cosmic time t, and the prime a derivative w.r.t. 
the scalar field 93. Eq. (10) then implies 

(fi^ = -2H'(p, (21) 

hence either 99 = or 

if = -2H', (p = m"H'. (22) 
After this is substituted back into the equations of motion we find |TT1 [TH [T3] 

V = 2,H^-2H'^, V' = 2H' {3H - 2H") . (23) 

It follows, that if) and H' can never vanish when V < 0, and conversely a configuration in 
which (p = H' = can be reached only in regions where V > 0. It also follows, that stationary 
points of V impose a particular constraint on -^[92]: 

V' = =^ H' = or 2H" = 3H. (24) 

In contrast, a stationary point of H will occur at a stationary point of V, unless H" is singular 
there and H'H" is finite and non-zero. In the latter case, eqs. (22) imply that = 0, (/3 7^ 0; 
hence if reaches an extremum and its trajectory in the (t, 99)-plane exhibits a turning point 
there. 



3 Explicit examples of scalar cosmology 

There are several ways to construct solutions for the equations of scalar cosmology, depending 
on the problem to be addressed. For example, one can search for models and initial conditions 
allowing for a specific type of solution, or one can try to find all solutions allowed by a given 
scalar potential and initial conditions. In this section we discuss examples of the first type; 
in the later sections we construct more generic solutions for specific potentials. 

Various cosmological scenarios, such as purely matter or radiation dominated universes, 
are described by a simple power law: 

a(t)=f-V, (25) 
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where r is some fixed reference time. Such types of behaviour can also be achieved in scalar 
field cosmology [H]; it follows from (25) that 



H 



a 
a 



n 
i' 



and the second equation ( 10 ) becomes 



-2H 



2n 



'2n 



t 



(26) 



(27) 



Clearly a power-law solution (25) is possible in the context of a regular scalar field model 



only for n > 0, i.e. non-contracting universes. Introducing a constant of integration r, the 

(28) 



solution of eq. (27) is O IH] 



ip{t) = V3(r) ± V2nln-. 

r 

Substitution of these results into the third eq. ^ now leads to 



V = 3H^ - -<j? 

2 ^ 



n{3n — 1) 



with Vo fixed by the requirement 



ip{T) _ n (3n - 1) 



(29) 



(30) 



Observe, that for < n < 1/3 the potential is negative definite, whereas for n > 1/3 it is 
positive definite. We conclude, that both positive and negative exponential potentials can 



allow for power-law solutions of the type ( 25 ) , but only in specific domains of non- negative 
powers n. We will show later, that flat contracting universes can arise in scalar cosmology for 



other types of potentials, although not with a simple power-law (25) for the scale factor. A 



different question to be addressed later is, what other solutions exist for exponential potentials. 

It goes without saying, that a constant Hubble parameter Ho results from constant ip = (pQ 
at an extremum of the potential: 



9? 



0, 



V'iipo 



0. 



(31) 



In scalar cosmology this necessarily represents a positive cosmological constant V 
Of course, a solution (31), as well as the Minkowski limit Hq = 0, can also arise as the final 



3Hl 



stationary state in a scenario in which the scalar field evolves dynamically from a higher value 
to end up at an extremum of the potential. 

Instead of specifying a certain evolution of the scale factor, one can also start from a spec- 
ification of the time-dependence of the scalar field. As an example, we construct a solution 
with an oscillating scalar field 

ip{t) = ipQCOsmt. (32) 

The rate of change of the field is 



The first eq. ( 22 ) then becomes 



-uiipo sin ojt = —uj \l ip^ — ip"^ . 



(33) 
(34) 
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with the solution 

^^l „„„„„„ , 1 , /,.2 ,^2 



H = Hq — arccos h - i^V'V ~ ■ ("^^) 



Here f^o is the initial value of H when tp = ipQ] indeed, substitution of (32) gives the time 
dependence of H as 

H{t) = HQ-\iplojH+\iplojsm.2ujt. (36) 
4 8 

The corresponding solution for the scale factor reads 

a{t) = aoe^«*-i '^l^^t^+h^lM"''^^^^^) _ (37) 

This represents a universe growing from very small size at large negative times, to a finite 
size around t = AHq/ {ojipfjf' , when it shows some oscillating behaviour, to contract again 
to arbitrarily small size for very large positive times. Finally, we can compute the scalar 
potential from which such behaviour follows: 



V = 3H^ - 2H' 



/2 



3 IHo — arccos — + ^ ^'Py'Po - f j - ^VPq-'P ) ■ 



(38) 



It is rather remarkable that such a complicated scalar potential can give rise to a simple 
periodic solution for the field (/^(t), and allows a complete solution for the scale factor. More 
importantly, we observe that this solution of the scalar cosmology equations clearly shows 
essentially reversible behaviour, and illustrates explicitly that the back reaction of the space- 
time curvature on the scalar field in the Klein-Gordon equation can not be interpreted off-hand 
as a dissipative friction term. 



4 More on exact solutions 

So far we have constructed potentials starting from a prescribed time-dependence of the scale 
factor a{t) or the scalar field ^{t). There is yet another way of finding solutions for scalar 
cosmology models, allowing to construct particular solutions as well as generic ones. This 
method starts not from a prescribed time behaviour of the cosmological degrees of freedom, 



but from postulating the relation H[lp\ in eq. (20). In this section we only consider particular 



cases with simple analytic solutions. The construction of generic solutions is discussed later. 
We start with the simplest non-trivial example, in which H is linear in the field ip: 

H = hQ + hnp, H' = hi. (39) 

It follows directly, that 

V = 3{ho + hipf -2hl (40) 
To diagonalize the mass term, we need to make the shift 

V' = (/? + ^ y = -2/if + 3/i?V^ H = hiil). (41) 

hi 
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Now hi is directly porportional to the mass: 



m^ = 6hj 1/ = + ^ m^^^ H = ^ip. (42) 

3 2 ^/Q 

The corresponding solutions for the field and scale factor are easily found: 



2 

m 



(43) 



This is not the only solution for the quadratic potential V in ( 42 ) , in fact it is a very special 
one: the only solution which exists for all times; but it pays to consider it a bit more in detail. 
A first observation is, that the rate of change of the scalar field ip is constant over the whole 
time domain {—oo, +oo). Thus there is no dissipation of kinetic energy, in spite of the fact 
that H does not vanish except at t = to- Secondly we observe, that the Hubble parameter 
H{t) is negative for times t > to; indeed this universe expands only during the epoch t < to 
from arbitrarily small scales to a maximal size when a{t) = oq for t = to, and contracts again 
to vanishingly small size at large positive time t > to- 

Actually, this behaviour for large times is generic for potentials with a negative minimum: 



Knm < 0. This follows from two general observations. First, eq. (10): 

implies that for regular kinetic terms of the scalar field the Hubble parameter is a non- 
increasing function of time, and is constant only at stationary points of the field evolution: 



if = 0. Second, as eq. ( 23 ) shows, any point where the potential is negative must satisfy 

V <0 2H'^ > 3H'^ > 0, (44) 

and therefore (p = —2H' can never vanish in the range where the potential is negative. As a 
result a negative minimum of the potential V, even if it is the absolute minimum, can never 
represent a stationary point of the dynamics; we conclude that the field never comes to rest 
at a negative value of the potential, and the Hubble parameter is a monotonically decreasing 
function of time as long as V is negative. This conclusion is in agreement with the general 
result of ref. [I6j . 

It is not difficult to construct more examples of exact solutions for polynomial potentials 
by a similar procedure. For example, taking 

H = ho + h2ip^, H' = 2h2ip, (45) 

we get a quartic potential 

V = Vo + '^ip' + ^^\ (46) 

with 

Vo = 3hl, = 12/10/12 - 16/i|, A = Uh^, (47) 
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which imphes 



4A / 3m2 ^ ^ 



^0 = ^1 + 



9 V 4A 

The corresponding particular time-dependent solutions of the Friedmann and Klein-Gordon 
equations are 

ifit) = ^(0)e-^*, a{t) = a(0)e\^*+i^'(°)(i-^"'"'), (48) 
with = 4A/3. The cosmology of this model was discussed in detail in ref. |13j . 

5 On a square root of the Friedmann equation 

In some cases it is possible to use the above procedure to construct the complete set of 
solutions in closed form. To simplify the discussion it is convenient to rescale the scalar field 
and define 

u{t) = (49) 



With this change of variable, the Friedmann equation (23) and the first field equation ( |22[ ) 
become 

H^-Hl = ^V, u = -3Hu, (50) 

employing the notation = dH/du. Now suppose the potential is positive definite: y > 0, 
in some domain of values of u. We can then introduce a function K[u) defined by 

H = ±^coshK, Hu = ±^{KuSmhK+'^coshK^. (51) 

Observe, that H can have either sign, but once the sign is fixed it cannot change anymore 
during the subsequent evolution of the universe; in the following we focus on positive H so 
as to describe an expanding universe. 



After substitution into eq. (50) and taking a square root, the function K is then seen to 
satisfy 

Ku + ^. cotanhK = ±1. (52) 

Observe that the equation is odd in K, hence in this equation the two sign choices are related 
by K — )• —K. Therefore it is sufficient to consider only the case with +1 on the right-hand 



side. An example is provided by the exponential potential (29) [T5]-[22]: 



leading to a very simple equation for K: 

+ ^ cotanh = 1 . (53) 

First consider the special cas^A = 2, for which 

Ku = l- cotanh K. (54) 



^The case A — —2 is obtained directly by the transformation u — >■ —u and K — > —K. 
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Observe, that cannot vanish anywhere, except in the Hmit K ^ ±00. The general solution 
of these equations is given implicitly by 



2K - e 



2K 



4(m - uq), 



(55) 



for some constant of integration uq . The Hubble parameter for an expanding universe is then 
determined by 

H = ly^SVo^ie^'' + 1) e-i(2^+-^"). 



The explicit time dependence can be obtained from the relation 



The two equations (56) and (57) can be combined to write 

1 



3H = -^e 



2K 



+ l)K = K 

It follows that there is a direct relation between K, u and a: 



u. 



a^e" ^ = constant. 



The constant defines a reference scale ao such that 



(56) 
(57) 

(58) 
(59) 

(60) 



Using this result one can eliminate K in terms of a and u. In addition, it also allows us to 
calculate the total expansion factor in some period of evolution, as expressed by the number 
of e- folds: 

N = \n— = \{K2- Ki-U2 + ui) = ^ (2K2 - 2Ki + e^^^ - e^^^) . 
Similar results can be derived for A = — 2. 

Having disposed of the cases for which = 4, we next turn to the generic case A^ 7^ 4. In 
terms of the initial condition Kq = K{uq) such that 



-2Kq/X 



the full solution is then given by 



1 + 



Ko 



1 



K+-ln 
2 



1 + 



-K 



Equivalently, 



1 + 



-K 



The corresponding expression for the Hubble parameter is 

,2K/X 



4 



{u - no) 



K 



6 



1 + 



K 



1-4 e 



(61) 

(62) 
(63) 

(64) 
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The pair of equations (62) and (64) represent the parametrized general solutions for {u,H), 
with time eliminated in favor of the parameter K. A well-known special solution of this kind 
is one for which K is constant: 



cotanh K 



A' 



Ku = 0, 



which requires |A| < 2. The Hubble parameter is then given by 

eA«/2 



H 



3 v^l-AV4' 



(65) 



(66) 



with opposite signs for an expanding or contracting universe. It is easy to check directly by 
substitution that this is a solution of the Friedmann equation. To construct the dynamics 
explicitly, observe that eq. ( 52 ) implies that 

fv 

Hu = \ — sinhi^, 



u 



-3Hu 



For the special solution ( 66 ) this leads to the results 



SAVo 



H 



^SVsmhK. 



1 



3A2 t-to' 



(67) 



(68) 



In the general case, by using eq. (62) one finds 
2k 



1 + 



K 



K 



(69) 



K 



This result can be used again to derive a direct relation between a, u and K. Indeed, eqs. 



(64) en (69) together imply 
3H = 



[l-l)e^ A 



K -ii 



(70) 



For A — )• ±2 this reproduces the results ^58h. The relation (59) now generalizes to 



constant 



3A/2 



(71) 



It is straightforward to extend the construction above to potentials which are negative definite: 
V < 0. This allows us to parametrize H as 



H 



sinhQ, Hu 



Vu 

Qu cosh Q + 1^ sinh Q 



(72) 



for some function Q{u). In contrast to the previous case, eq. (51), here there is no need of 
a sign choice, as it can be absorbed in the sign of Q. Moreover, H can change sign during 
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the evolution of the universe, in case Q switches sign. We have noted before, that this is a 
fundamental difference between strictly non-negative potentials and potentials taking negative 
values in some domain of scalar field values. 

By taking a square root, the Friedmann equation becomes 



Vu 

Qu + ^ tanhQ = ±1, 



(73) 



where again the two sign choices are related hy Q ^ —Q and we can restrict ourselves to the 



positive sign without loss of generality. Using the example of the exponential potential (29), 
with negative amplitude: 



V = Voe 



Xu 



Vo<0, 



eq. (73) becomes 



For 7^ 4 the solution is 



Q„ + - tanhQ = 1- 



l-^]{u-uo) 



or equivalently 



= (x-t)("-"o) ^ p2Q/X 



1 + 



Q 



+ 1 



(74) 



(75) 



(76) 



In this way we again construct a parametrized solution for the pair {u,H). For |A| > 2 there 
exists another special simple solution, with constant Q: 



tanhQ = 

The corresponding Hubble parameter is 

H = ±^ 



Vo 
3 



Qu = 0- 



A2 



(77) 



(78) 



For the general solution 



H 



sinhQ, n = — Y^|3y[cosh Q, 



(79) 



and from (74): 



Q = uQu = \/|3y| sinhQ - coshQ 



3A 



H + ii. 



(80) 



This implies a relation between Q, u and a similar to (71): 



\ 3A/2 



(81) 
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6 Series expansions: the regular case 



For general potentials, even if one cannot produce exact solutions, one can always construct 
solutions for scalar cosmology based on a series expansion method. In this section we consider 
the regular case, in which ip(t) is a single- valued function of time in some finite time domain. 
We have already observed before, that in single-scalar cosmology H{t) is a non-increasing 
function of time, and therefore it can be represented by a well-behaved function H[(p{t)] in 
the time domain considered. Later we will also consider solutions in a time domain in which 
ip{t) has a turning point, and the equation for H[ip] becomes double-valued. 

As in sect.jHj it is convenient to work with a rescaled scalar field u = y^3/2 and a Hubble 
parameter -ff[n], satisfying the equations (50). Let uq be a point in the regular domain; then 
we can develop H in a power series 



H = hn{u - Up)" = hp + hi{u - Up) + h2{u - u^f + 

n>0 

Hu = y^(n + l)hn+i{u - uqY = hi + h2{u - up) + ... 



(82) 



n>0 



Clearly, in this case 



Vu = 6Hu{H-Huu) = <^ H^ = 0. (83) 

i.e., H can have a stationary point only at an extremum of the potential; in all other points 
H[u] is a strictly monotonically decreasing function of time, hence a monotonic function of 
u, decreasing or increasing for positive or negative slope of u{t), respectively: 

Hu<0, 

Vu^O { (84) 

Hu>0. 




For example, a quadratic potential 



V = e + —u', (85) 



has a single minimum at u = 0. Hence H can have a stationary point Hu = only at u = 0, 
and 

Hu = <^ ^' = |- (86) 

This condition can only be fulfilled if V{0) = e > 0, and the stationary solution is a Minkowski 
space for 

n = 0, e = H = 0, (87) 
whilst it becomes a de Sitter space if 

u = 0, e = 3H'^ > 0. 



There is no stationary solution for e < 0. Note, that all other solutions for any e are non- 
stationary, with u = —3Hu 7^ 0. 

Returning to the series expansion (82), and a similar expansion for the potential: 

V = ^w„(u - no)" = vo + vi{u - Up) + V2{u - up)"^ + (89) 

n>0 
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eqs. (50) give rise to the infinite set of relations 

n 

^ [hkK-k -{k+l){n-k + l)/ife+i/in-fc+i] = y , (90) 

A:=0 

for all non- negative integers n = 0,l,2,..., and to 

u = -2,Y,{ri + ^)hn+i{u-UQT. (91) 

n>0 



(92) 



The first few relations (90) in explicit form are 

n = : VQ = 2>hl- 3hj, 

n = 1 : vi = 6/ii {ho — 2/12) , 

n = 2 : V2 = 3hi {hi - 6/13) + 6/12 {ho - 2/12) , 

n = 3 : V3 = 6hi (/i2 — 4/14) + 6/13 (/iq — 6/12) • 

It follows that either hi = Hu{uo) = 0, which can happen only if vi = Vu{uo) = and 
u{uo) = 0, or /ii / and 

h2 u2 , ^0 r ho ( ^ Vi 

hQ = hi + —, h2 



3 ' 2 y (jhohi 

^ /to fs /to / \ ( ^0^1 ^1 9 

^~ 24 24/11 72/i2 V 4/io/iiy V ^1 S/i? 

Using these last results, the equation for the scalar field becomes 

Thus the complete solution is given in terms of two parameters ho = H{uo) and hi = Hu{uo), 
representing the initial conditions of the cosmology. 

A solution with hi = exists only if = and uq is an extremum of V; then eqs. (92) 
reduce to 

vo = 3/iq, Vl = 0, 

V2 = 6/io/i2 — 12/i|, V3 = dhohs — 36/12/13, (95) 

= 6/io/i4 - 48/i2/i4 + 3/ti - 27/i|, ... 

Therefore, if ^2 7^ we have Huu = 2/i2 7^ 0, and the point mq is a point of inflection of 
where the field comes to rest (either momentarily or permanently). Obviously such solutions 
are very special, if they exist at all for some given potential. It requires a trajectory H[u] 
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to reach a local extremum of the potential at zero velocity. An example is provided by the 
special solution of the quartic potential in eq. (45) and following, which comes to rest at the 
minimum 999 = of the potential ( 46 ) . 

From the results above one can also estimate the total expansion factor of the universe between 
two times {ti,t2), as given by the number of e-folds. The central result is, that 



N 



Hdt 



H 

3Hu 



-du 



du ^0 + hi{u - uo) + h2{u - uo)'^ + .. 
hi + 2/i2(u - no) + 3/i3(n - no)^ + 



For the generic case hi the result is again a power series expansion 

2 



k>l 

with coefficients 



h( 



u 



uo) + n2{u-uo) +n3(n 



uof + 



ni 



Iho 

'3 hi 



n2 



1 1 hQh2 
'6^3 h\ ' 



(96) 



(97) 



1 /i2 _^ 1 /iq/is 



4 /lo/io 



^ hi 3 hi 9 /if ' 



(98) 



n4 



1 /is Ihoju 
6/11 3 hi 



Qhi 



hoh2h3 



hi 

For the special case hi = one gets in stead an expansion 



+ 



2 hohl 



3 h\ 



with 



= [no ln(n - no) + ni(n - no) + n2(n - no)^ + n3(n 

1 /lo/l3 



uof + ..]]. 



no 



Iho 
'6/12' 



ni 



4 hi ' 



(99) 



n2 



1 hphj 
12 6 hi 



3 hoh 



16 hi 



(100) 



ns 



36 h2 



5 /lo/i, 



0"5 



36 hi 



1 /lo/i3/i4 3 /iq/Is 

3 /i3 16 hi 



7 Applications 

In this section we apply the general results above to simple quadratic potentials (85). The 
simplest models are those with e = 0, which have a Minkowski minimum H = at u = 0. It is 
most convenient to expand around the minimum no = 0. Then there is only one non- vanishing 
term in the potential 



771 

V2 = ^, Vn = 0, n = 0,1,3,... (101) 
As a result we get for the non-stationary solutions which all have hi 7^ 0: 



ho 
hi 



±1, 



h2 
hi 



1 



h3 
hi 



1 



1 



m 
9/if 



/14 

h'l 



1 

"24 



1 + 



2m^ 
9/1? 



(102) 
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The power series expansion for H[u] then takes the form 

,2X 1 



H 



±h 



1 5 1 



2 6 
±hi [e±" + 0(mV/i?)] 



1 



m 



24 



1 + 



2m^ 
9J4 



(103) 



This result was to be expected, as in the hmit — )• the potential vanishes and the solutions 
of the Friedmann equation ( 50 ) become pure exponentials. For the equation of motion of the 
scalar field we get similarly 



u 

3h~i 



l±u + 



1 



1 



m 
9)^ 



u 



+ ... = + O {my hi) . 



(104) 



u 



Fig. 1: Physical domain in the («, _f/)-plane for e = as determined by eq. (|105| 



The domain of validity of this series expansion is restricted by the requirement that H[u] is 
single valued. Now u{t) can have a turning point only where u = —3Hu = 0. Therefore the 
locus of potential turning points is 



V 
3" 



mu \ 2 



(105) 



for the case at hand. Thus the (u, H)-p\ane is divided in four sectors by straight lines solving 



(103): 



eq. ( 105 ), and only the upper and lower quadrants in fig. 1 are allowed regions for the solutions 

(106) 



mu 


< H[u] < 


mu 






3 




3 



There is no solution crossing from the upper to the lower quadrant. Indeed, at the point 
u = H = where the lines cross, ii = = and any solution passing through this point 
must have ho = hi = 0, hence Huu{0) = 2/12 = 0. It follows that no solution can pass 
from positive to negative H, and the two sets of solutions are strictly separated. The only 
exceptional case is the Minkowski solution represented by the point at the origin. 
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u 



Fig. 2; Physical domain in the (u, i?)-plane for e > as determined by eq. (|109 



Next we turn to the case e > 0, which has a stationary solution u 



representing 



de Sitter space, with a cosmological constant given by eq. (86). The solution for H is now 
slightly modified to 



ho 
hi 

h 

hi 



with the result 



H 



±. 1 + 



£ 


h2 


_ ho 


3/if' 


hi 


~ 2hi 


Tin? \ 




1 










h~i 


= =^24 



±hi 



1 + 



4f^i- 



e 
3)^ 



1 + 



2m^ 



e 9 1 



m 



u 



(107) 



+ 



1 + 



1 + 



2m' 

2 

m 
9^ 



(108) 



In the {u, -ff)-plane the domain of validity of these approximations is restricted by the branches 
of the hyperbola 

" -{-) -I- (™' 

shown in fig. 2. On this hyperbola u = Hu = 0. The two domains of allowed positive and 
negative if- values are separated by a gap of size = 2y^e/3. No cross-over is possible, and 
the solutions describe only permanently expanding or permanently contracting universes. 
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u 



Fig. 3: Physical domain in the (u, //)-plane for e < as determined by eq. ( |112[ ). 

Finally, we consider the case e < 0. If there would be a solution u = u = 0, this would 
give rise to an anti-de-Sitter space. However, no such solution exists: in the domain V < 
solutions of the Klein-Gordon equation are always dynamical: 

u = -3Hu, Hl = -- >0, for all V <0. (110) 

3 



In particular 

hence it is guaranteed that hi ^ 0. The solution for H and for ii is formally the same 



^? = /io + y>0, (111) 



as in eqs. (109), except that on has to replace e = — The restriction imposed by the 



single- valuedness of H[u] now becomes 

which is a hyperbola with branches in the left and right quadrants, leaving an opening on the 
n-axis in the interval 



el \/\e\ 
<u < 



m m 

This hyperbola is shown in fig. 3; passing through the allowed interval on the u-axis, H 
can cross from positive to negative values. As H{t) is a non-increasing function, this will 
eventually happen and an expanding universe will turn into a contracting universe. We have 



seen this behaviour already in the example given by eqs. (42), (43). We return to this point 
for a fuller discussion in the next section. 



8 Series expansions: turning points 



As eq. (10) shows, in single scalar cosmology H < 0. Now the scalar field u 

Hu 



have turning points where n 
notation 



0, but il 



K = (H - H.. 



Vu 7^ 0. Adapting eq. (23) to the present 

(113) 
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hence such a turning point occurs if Huu is singular in such a way that at this point 

< \HuHuu\ < oo. (114) 
As discussed in the previous sections, in the {u, //)-plane turning points he on the curves 

V{u) = 3H^{u), 

which bound the domain of physicaUy ahowed values. In fact all points on these boundary 
curves are turning points, unless Vu = 0, i.e. at a local extremum of the potential. In the latter 
case a solution with ii = = can exist only if this extremum occurs at a non-negative 
value of V. 

In the neighborhood of a turning point u{t) takes identical values before and after the 
turning point; but as -u = —3Hu ^ away from the turning point, and therefore H < both 
before and after, H[u{t)] necessarily becomes double-valued there. As a result, the power 
series expansion studied in sect. |6]can not be used in this neighborhood. 

This double- valuedness can be resolved by a reparametrization of the field. For definiteness 
it is convenient to consider a point where u reaches a maximum Um', then a new dynamical 
variable r](t) can be introduced such that in a sufficiently small but finite time interval around 
the turning point 

u{t) = Um-v\t)- (115) 

At the turning point r] = 0, and the evolution of u{t) can be parametrized by a monotonically 
increasing function ?y(t), the negative and positive values of rj corresponding to u before and 
after the turning point, respectively. In case Um were a minimum, we could similarly define 



u{t) = Um + rf{t). 



(116) 



However, for our discussion we will assume a maximum and use eq. (115). 

We first re-express our dynamical equations in terms of the new field r]{t). Using the 
definition (115) and the notation Hn = dH/drj, it is straightforward to derive the equations 



f7-2 1 fr2 



V[u{t])], f] 



These equations can be used to develop a new power series expansion 



H 



n>0 



n>0 



[111) 



(118) 



Now expressed in rj the condition (|114|) at the turning point, where Hu = 0, reads 

1 



V\ri=0 



0, 



< 



3 Hj^Hr/ri 



< OO. 



ri=0 



(119) 



In terms of the expansions (118) this is translated as 

gi= 92 = 0. 

As a result 



(120) 



H = go + g^rf + g^rf + bg^rf..., 



Hr, = Sgsi]^ + Ag^rf^ + bg^v'^ + 



(121) 
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Also, assuming that the potential V has a power series expansion (89) with uq = Um, the 
potential has the expansion 

V = ^{-irvnV^''- (122) 



n>0 



Then eqs. (117) imply for the coefficients gn- 



2 _ ^0 

50 -y, 



9a 



90 



91= 92 = 

93 



9l 



95 



{2vo - 9v2) 



(123) 



96 



290 
405^1 



-2vo + 9v2) , 97 



253 
21 



1 + 



9^3 



+ 



1 



(2^0 - 9v2) i2vo + 15v2) 



etc. Observe, that these equations require t'o ^ and vi > 0. The first condition implies that 
turning points only occur at non- negative values of the potential V. The second condition is 
a direct consequence of our choice to consider a turning point which is a maximum of u: 

Vu{Um) = f 1 > 0. 

The equations for H[rj] and r]{t) now become 



H 



2 /^Tf 3 1 4 9V2 - 2vo 5 

1 \ — 1] -\ — rj -\ 7] + ... 

3 V ^^0 3 Ab^/vQVi 



(124) 



3f 1 



3 V ^1 



2vo — 9v2 
18vi 



rf + .. 



Note, that in order to get a monotonically increasing ?7(t), we have to take the negative square 
root for (73: 

2 

In the present formulation one can derive yet another formula for the total expansion 
factor in a given time interval (ti,t2): 



Hdt 



drj 



^ , 90 + 93V^ + 9'iV^ + 95V^ + 
arj ■ 



(125) 



1 353 +45-47/ + 555 7?2 + ... 

The result can be written in a series expansion as 

i2 



N = [vilf] + U2lf + 1^3 rf + 1^41]^ + ....] ^ , 



(126) 



with coefficients given by 

450 
'953' 



1^2 



20 9095 64 gogi 

^3 ~ ^ — 2 ^ — 3^' ^4 

81 g^ 243 



8 goff4 
27 gj ' 



12 goge, _ 40 ffog4g5 64 gosi 



(127) 



9 9 gi 81 53^ 



243 
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Note again, that by taking the negative value for the first coefficient becomes positive: 




(128) 



9 Quadratic potentials 



The general description of turning points developed in sect. [8] can be illustrated with the 
example of quadratic potentials considered earlier 



m 



2 



We consider a solution u{t) which rolls down the potential from negative values to positive 
values, reaching a turning point at some positive Um > 0, where Um = and Vu{um) > 0. 



Using the parametrization (115) the potential is expressed as 



2 4 2 
V = vo-vir] +V2r] , vq = e + —u^, vi = —^Um, V2 = —. (129) 

The requirement > is satisfied automatically, whilst the condition > is non-trivial 
only if e < 0, requiring Um to be in the domain of V{um) > 0. As before we distinguish the 
cases e > and e < 0. 



For non-negative e the coefficients gn in eq. (123) become 




2m r- — 

2 2' 9z = --prV^Um, 

m^n£, 9 



54="^Jl + 4V' g5= "^ fl-'^-^V (130) 




2ul 2e \ 
9 3m2 J 



The singularity of g^ for Um — ^ indicates, that the only consistent solution with n = at 
the point « = is the one with constant H = go = y^e/3, corresponding to a de Sitter space 
for e > 0, or a Minkowski space for e = 0. This is a stationary solution, rather than a turning 
point. 



For negative e turning points can occur, and eqs. (130) still hold, if > 3|e|/m2. 
Therefore, if the scalar field starts at a large enough value in the region V{u) > 0, it can roll 
down the potential and oscillate, meeting a number of turning points, until it can no longer 
escape from the region V < and the universe starts to contract back to infinitely small size. 
The point is, that after each turning point the Hubble parameter will continue to decrease: 
H < 0, until it finally crosses over into the region of negative H and the contraction phase 
starts. 
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10 Discussion and conclusions 



The main concern of this paper is the cosmological evolution of spatially flat universes driven 
by a single scalar field. Such a scenario may be relevant for the present epoch in the evolution 
of our accelerating universe, and it may have a bearing on the vey early universe going 
through an epoch of inflation. Observations of the CMB require an inflationary expansion by 
at least 60-70 e-folds, which can happen only if the universe spends a relatively long time in 
a phase with large Hubble parameter. This requirement is usually expressed by the slow-roll 
condition; therefore it is of interest to study this condition in the present context of single 
scalar-field cosmology. 

The accelaration of the universe can be expressed as 



a 



H + H"^ = (1- e) , (131) 
a 

where the slow-roll parameter is defined by 

^=^- (132) 

Thus an accelerated expansion requires < e < 1. Now we can combine the two inequalities 

> 3Hl, 3H^ -V = ml > 0, (133) 

to translate the slow-roll condition to a double bound on H^: 

V o V 

Clearly these bounds can be satisfied only in a region where V > Also, the bound is always 



satisfied at and near a turning point where = 0. In terms of our series expansion (121) 
this is at = 0, where = y/3. In order to estimate the total expansion factor in the 
slow-roll domain (134) near a turning point, we also ought to find a value for r/ at = V/2 



for a given specific potential. In general however, we can use the observation that in Planck 
units Um and rj must satisfy 

if < u.m < 1, (135) 

and as a first approximation 

9 93 Jo V ' y 3V^^i ^ ^' 

where r]m is the upper limit for r] where H = V/2. The condition on the number of e-folds 
for inflation then becomes 

^ ~ 10^ 

Vl 

However, significant modification of this estimate may result for specific potentials j23j . 

The methods we have used in this paper to find solutions for the the cosmological equations 
of a spatially flat universe driven by scalar fields relies heavily on the fact that we have assumed 
a single field to drive the cosmological expansion. This allows us to replace time by the field 
as the evolution parameter. Actually eq. ( |l2[ ) suggests another option, taking X'^ = \/61na as 
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evolution parameter; this might be more readily generalizable to the case of many scalar fields 
|24t [9]. At present such a modification of the methods presented here is under investigation. 
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